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The fact that the simplest modern cosmological theory, standard Cold Dark Matter (sCDM), almost fits all 
available data has encouraged the search for variants of CDM that can do better. Here we discuss what are 
perhaps the two most popular variants of CDM that might agree with the data: ACDM and CHDM. While the 
predictions of COBE-normalized ACDM and CHDM both agree well with the available data on scales of ~ 10 to 
100 h~ l Mpc, each has potential virtues and defects. ACDM with fio ~ 0.3 has the possible virtue of allowing 
a higher expansion rate Ho for a given cosmic age to, but the defect of predicting too much fluctuation power 
on small scales. CHDM has less power on small scales, so its predictions appear to be in good agreement with 
data on the galaxy distribution, but it remains to be seen whether it predicts early enough galaxy formation to 
be compatible with the latest high-redshift data. Also, several sorts of data suggest that neutrinos have nonzero 
mass. And two very recent observational results favor high cosmic density, and thus favor Q = 1 models such 
as CHDM over ACDM — (1) the positive deceleration parameter go > measured using high-redshift Type la 
supernovae, and (2) the low primordial deuterium/hydrogen ratio measured in two different quasar absorption 
spectra. If confirmed, (1) rules out a cosmological constant large enough to help significantly with the Ho-to 
problem; while (2) suggests that the baryonic cosmological density is at the upper end of the range allowed by 
Big Bang Nucleosynthesis, perhaps high enough to resolve the "cluster baryon crisis" for Q = 1 models. We try 
to identify "best" variants of both ACDM and CHDM, and discuss critical observational tests for both models. 



1. INTRODUCTION AND SUMMARY 

"Standard" ft = 1 Cold Dark Matter (sCDM) 
with a near-Zel'dovich spectrum of primordial 
fluctuations [1] seemed to many theorists to be 
the most attractive of all modern cosmological 
models. But while sCDM normalized to the am- 
plitude of cosmic microwave background fluctu- 
ations observed by COBE is a good fit to most 
large-scale observations, it is inconsistent with 
data on smaller scales [2]. Here we discuss what 
are perhaps the two most popular variants of 
sCDM that might agree with the data: ACDM 
and CHDM. The linear matter power spectra for 
these two models are compared in Figure 1 with 
the real-space galaxy power spectrum obtained 
from the two-dimensional APM galaxy power 
spectrum [3]. The ACDM and CHDM models 
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essentially bracket the range of power spectra in 
currently popular cosmological models which are 
variants of CDM. 

ACDM cosmological models with a positive 
cosmological constant (A > 0) and Slo = 1 — S^a ~ 
0.3, where £l\ = A/(3Hq), have been advocated 
[4] because they allow a larger Hubble constant 
Ho for a given age of the universe to, and they 
predict a larger fraction of baryons in galaxy 
clusters, than SI = 1 models. Early galaxy for- 
mation also is considered to be a desirable fea- 
ture of these models. But early galaxy forma- 
tion implies that fluctuations on scales of a few 
Megaparsecs spent more time in the nonlinear 
regime, as compared with Cold + Hot Dark Mat- 
ter (CHDM) models. As has been known for 
a long time, this results in excessive clustering 
on small scales. We have found that a typical 
ACDM model with H = 70 km s" 1 Mpc" 1 , 
Sl = 0.3, and cosmological constant A such that 
fi A = A/(3if 2 ) = normalized to COBE on 

large scales (this fixes erg « 1.1 for this model), 
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is compatible with the number-density of galaxy 
clusters, but predicts a power spectrum of galaxy 
clustering in real space which is much too high 
for wavenumbers k = (0.4 — l)/i/Mpc [5]. This 
conclusion holds if we assume either that galaxies 
trace the dark matter, or just that a region with 
higher density produces more galaxies than a re- 
gion with lower density. One can see immediately 
from Figure 1 that there will be a problem with 
this ACDM model, since the APM power spec- 
trum is approximately equal to the linear power 
spectrum at wavenumber k « 0.6/i Mpc -1 , so 
there is no room for the extra power that nonlin- 
ear evolution will certainly produce on this scale 
(see Figure 1 of Ref. [5] and further discussion 
below). The only way to reconcile the model 
with the observed power spectrum is to assume 
that some mechanism causes strong anti-biasing 
— i.e., that regions with high dark matter den- 
sity produce fewer galaxies than regions with low 
density. While theoretically possible, this seems 
very unlikely; biasing rather than anti-biasing is 
expected on small scales [6]. 
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Figure 1. Power spectrum of dark matter for 
ACDM and CHDM models considered in this pa- 
per, both normalized to COBE, compared to the 
APM galaxy real-space power spectrum. 



Our motivation to investigate this particular 
ACDM model was to have Ho as large as might 
possibly be allowed in this class of models, which 
in turn forces Slo to be rather small in order to 
have to > 13 Gyr. There is little room to lower 
the normalization of this ACDM model by tilting 
the primordial power spectrum P p (k) = Ak" p , 
i.e., assuming n p significantly smaller than the 
"Zel'dovich" value n p = 1, since then the fit to 
data on intermediate scales will be unacceptable 
— e.g., the number density of clusters will be too 
small [5]. Tilted ACDM models with higher Slo, 
and therefore lower Ho for to ^ 13 Gyr, appear to 
have a better hope of fitting the available data, 
based on comparing quasi-linear calculations to 
the data [5,7]. But all cosmological models with a 
cosmological constant A large enough to help sig- 
nificantly with the Ho-age problem are in trouble 
with observations providing strong upper limits 
on A [9]: gravitational lensing [8], number counts 
of elliptical galaxies [10], and especially the pre- 
liminary results from measurements of the cos- 
mological deceleration parameter go using high- 
redshift Type la supernovae [11]. 

CHDM cosmological models with SI = 1 mostly 
in cold dark matter but with a small admixture 
of hot dark matter (light neutrinos with a total 
mass of about 5 eV contributing Q, v « 0.2 for 
h = 0.5, where h = H o /(100 km s _1 Mpc" 1 )) are 
a good fit to much observational data [12,13] — 
for example, correlations of galaxies and clusters 
and direct measurements of the power spectrum 
P(k), velocities on small and large scales, and 
other statistics such as the Void Probability Func- 
tion (probability Po(r) of finding no bright galaxy 
in a randomly placed sphere of radius r) . We had 
earlier shown that CHDM with Q, v = 0.3 predicts 
a VPF larger than observations indicate [14], but 
new results based on our Q, v = 0.2 simulations 
in which the neutrino mass is shared equally be- 
tween two neutrino species [12] show that the 
VPF for this model is in excellent agreement 
with observations. However, our simulations [5] 
of COBE-normalized ACDM with H = 70 and 
Sl = 0.3 lead to a VPS which is not compatible 
with the data [15]. 

Moreover, there is mounting astrophysical and 
laboratory data suggesting that neutrinos have 
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non-zero mass [12,16]. The analysis of the LSND 
data through 1995 [17] strengthens the signal 
for —?> v e ; comparison with exclusion plots 
from other experiments implies a lower limit 
Am 2 e = \m(i' fi ) 2 — m{v e ) 2 \ > 0.2 eV 2 , imply- 
ing in turn a lower limit m v > 0.45 eV, or 
Q. v £ 0.02(0.5//i) 2 . This implies that the con- 
tribution of hot dark matter to the cosmological 
density is larger than that of all the visible stars 
(£1* « 0.004 [18]). The Kamiokande data showing 
that the deficit of atmospheric muon neutrinos in- 
creases with zenith angle suggests that — > v T 
oscillations occur with an oscillation length com- 
parable to the height of the atmosphere, implying 
that Ara^ ~ 10" 2 eV 2 [19] — i.e., that and v T 
are nearly degenerate in mass, consistent with the 
hot dark matter mass being shared between two 
neutrino species. Because free streaming of the 
neutrinos damps small-scale fluctuations, even a 
little hot dark matter causes reduced fluctuation 
power on small scales and requires substantial 
cold dark matter to compensate; thus evidence 
for even 2 eV of neutrino mass favors large £1 and 
would be incompatible with Sl c as small as 0.3 
[12]. Allowing Slj, and the tilt to vary, CHDM 
can fit observations over a somewhat wider range 
of values of the Hubble parameter h than stan- 
dard or tilted CDM [20], especially if the neutrino 
mass is shared between two or three neutrino 
species [12,21]. However, as for all SI = 1 models, 
h larger than about 0.5 conflicts with to > 13 Gyr 
from globular cluster and white dwarf cooling age 
estimates. Another consequence of the reduced 
power on small scales is that structure formation 
is more recent in CHDM than in ACDM. This 
may conflict with observations of protogalaxies 
at high redshift, although the available evidence 
does not yet permit a clear decision on this (see 
below). The evidence from preliminary data of 
a fall-off of the amount of neutral hydrogen in 
damped Lyman a systems for z J> 3 [22] is in 
accord with predictions of CHDM [23]. 



2. CLUSTER BARYONS 

Primack has recently reviewed the astrophysi- 
cal data bearing on the values of the fundamental 



cosmological parameters, especially Slo [9]. One 
of the arguments against SI = 1 which seemed 
hardest to answer was the "cluster baryon cri- 
sis" [24]: recent X-ray observations of clusters 
show that the abundance of baryons, mostly in 
the form of gas (which typically amounts to sev- 
eral times the total mass of the cluster galaxies), 
is as much as 20% if h is as low as 0.5. For the 
Coma cluster the baryon fraction within the Abell 
radius is 

ft = 4r^> °- 009 + o.050/r 3 / 2 , (l) 

iVltot 

where the first term comes from the galaxies and 
the second from gas. If clusters are a fair sample 
of both baryons and dark matter, as they are ex- 
pected to be based on simulations, then this is 2-3 
times the amount of baryonic mass expected on 
the basis of BBN in an £1 = 1, /i « 0.5 universe, 
though it is just what one would expect in a uni- 
verse with Slo ~ 0.3. The fair sample hypothesis 
implies that 

A recent review of gas in a sample of clus- 
ters [25] finds that the baryon mass fraction 
within about 1 Mpc lies between 10 and 22% 
(for h = 0.5; the limits scale as h~ 3 l 2 ), and ar- 
gues that it is unlikely that (a) the gas could be 
clumped enough to lead to significant overesti- 
mates of the total gas mass — the main escape 
route considered in [24] (cf. also [26]). If SI = 1, 
the alternatives are then either (b) that clusters 
have more mass than virial estimates based on 
the cluster galaxy velocities or estimates based 
on hydrostatic equilibrium [27] of the gas at the 
measured X-ray temperature (which is surprising 
since they agree [28]), or (c) that the usual BBN 
estimate £16 k, 0.05(0.5//i) 2 is wrong. It is inter- 
esting that there are indications from weak lens- 
ing [29] that at least some clusters may actually 
have extended halos of dark matter — something 
that is expected to a greater extent if the dark 
matter is a mixture of cold and hot components, 
since the hot component clusters less than the 
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cold [30,31]. If so, the number density of clus- 
ters as a function of mass is higher than usually 
estimated, which has interesting cosmological im- 
plications (e.g., (78 is a little higher than usually 
estimated). It is of course possible that the solu- 
tion is some combination of alternatives (a), (b), 
and (c). If none of the alternatives is right, then 
the only conclusion left is that Slo ~ 0.33. The 
cluster baryon problem is clearly an issue that 
deserves very careful examination. 

It has recently been argued [32], CHDM models 
are compatible with the X-ray data within ob- 
servational uncertainties of both the BBN pre- 
dictions and X-ray data. Indeed, the rather 
high baryon fraction £lb « 0.1(0.5//i) 2 implied 
by recent measurements of low D/H in two high- 
redshift Lyman limit systems [33] helps resolve 
the cluster baryon crisis for £1 = 1 — it is es- 
cape route (c) above. With the higher Slo im- 
plied by the low D/H, there is now a "baryon 
cluster crisis", for low Slo models! Even with a 
baryon fraction at the high end of observations, 
fb ^ 0.2(/i/0.5) _3/ ' 2 , the fair sample hypothesis 
with this Sl 6 implies Sl ^ 0.5(/i/0.5) _1/2 . 

3. ACDM 

Although ACDM was considered from the be- 
ginning of CDM as a viable alternative to the 
standard SI = 1 version [1], the first N-body sim- 
ulations [34] revealed a problem with the model: 
the correlation function is too steep and too large 
at small scales. For example, Efstathiou et al. [4] 
found that the model with Slo = 0.20 nicely 
matches the APM angular correlation function 
w(6) on large angular scales (6 > 1°), but dis- 
agrees by a factor of 3 with APM results on 
smaller scales. The disagreement was not consid- 
ered serious because "the models neglect physics 
that is likely to be important on small scales 
where £(r) > 1". 

While it is true that one cannot reliably take 
into account all physics of galaxy formation, 
that does not mean that we can manipulate the 
"galaxies" in models without any restrictions. 
For example, one of the ways out of the prob- 
lem would be to assume that places with high 
dark matter density somehow are less efficient in 



producing galaxies. But Coles [35] showed that 
if the dark matter density has a Gaussian dis- 
tribution and the number density of galaxies is 
any function of local dark matter density ("lo- 
cal bias" ) , then the correlation function of galax- 
ies cannot be flatter then the correlation function 
of the dark matter. Coles' results are formally 
applicable only for a Gaussian density distribu- 
tion function, which is not the case for real mod- 
els at a nonlinear scale. But it illustrates our 
point: just introducing antibias does not neces- 
sarily save the situation. In the case of a Gaus- 
sian distribution it actually makes it worse. Our 
simulations [5] indicate that Coles' results might 
be true even for realistic non-Gaussian distribu- 
tions. This may help explain why we cannot sat- 
isfy APM w(6) constraints simultaneously at large 
and small scales with any local bias: because the 
correlation function of dark matter in the ACDM 
model was already too steep at small scales. (We 
can reduce small-scale w(6), but this will ruin it 
on large scales.) 

Because our numerical and Coles' analytical 
results were obtained for local bias, one might 
think that in order to reconcile the model with 
observations we need to appeal to nonlocal ef- 
fects [36] which might result from photoionization 
of the ISM by UV photons produced by quasars, 
AGNs, and young galaxies, or propagation of 
shock waves produced by multiple supernovae in 
active galaxies. Actually, neither effect is very ef- 
ficient in suppressing star formation in high den- 
sity environments. UV radiation heats gas only 
to few tens of thousands of degrees. This does af- 
fect the formation of small galaxies and delays the 
time of formation of the first stars in large galax- 
ies, but it cannot change a large galaxy [37,38]. 
If the gas falls into the gravitational potential of 
a normal-size galaxy with effective temperature 
of about 10 6 K, it does not matter much if it was 
ionized and preheated. Even if a strong shock is 
produced by an active galaxy, it is difficult to de- 
liver the shock to a nearby galaxy. Because galax- 
ies are formed in very inhomogeneous environ- 
ments, a shock produced by one galaxy will have 
a tendency to damp its energy into a local void, 
not to propagate into a dense area where another 
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galaxy is forming. Numerical hydro+N-body sim- 
ulations by Yepes et al. [39], which incorporate 
effects of UV radiation, star formation, and su- 
pernovae explosions, do not show any antibias of 
luminous matter relative to the dark matter. 

In our recent study of ACDM [5], we tried to 
avoid complicated questions about effects of star 
formation. We used several Particle Mesh dissipa- 
tionless simulations of various box sizes and reso- 
lutions to estimate the nonlinear power spectrum 
of dark matter, and reinforced the old conclusion 
that it is not compatible with the observed clus- 
tering of galaxies. Our results are shown as the 
solid line in Figure 2; as expected, the nonlinear 
spectrum has much extra power on small scales. 
Thus, this ACDM model, which has a formal bias 
parameter b ~ 1, must have extra (anti-) bias. 
Then we put lower limits on the possible power 
spectrum of galaxy clustering. This includes two 
steps: 

i) All regions with mass less than 5 x 
10 9 /i -1 Mq are assigned zero luminosity 
because each such region does not have 
enough mass to produce a galaxy luminous 
enough to make its way into the CfA or 
the APM catalogs. It also raises the power 
spectrum by a factor of four by removing 
the ~ 50% of mass in low-density regions. 

ii) We assume that in the remaining, higher 
density regions, the number density of 
galaxies does not depend on the density of 
the dark matter. 

These purposely unrealistic assumptions suppress 
the number of galaxies identified in groups and 
clusters, for example. Because galaxies in high 
density regions are expected to be more clus- 
tered than galaxies in the field, this scheme gives 
a lower limit on the galaxy clustering predicted 
by the model. It also gives a significant antib- 
ias for galaxies. But even with this antibias, 
the resulting lower limits on the power spectrum, 
which are shown in Figure 2 as the dashed curve 
(from our highest-resolution simulation) and dot- 
dash curve (from a larger box but lower resolu- 
tion simulation), are 2-3 times higher than CfA 
estimates [40] and 3-4 times higher than APM 



results [3]. (The CfA estimates were obtained 
by converting the redshift-space power spectrum 
to real space by comparing redshift-space and 
real-space power spectra from a large-box low- 
resolution ACDM simulation; both because of the 
uncertainties inherent in this procedure and be- 
cause the CfA survey is more shallow, we regard 
the APM results as being probably more reliable.) 
A more conventional galaxy identification method 
(e.g., high peaks) would imply even larger dis- 
crepancies. 
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Figure 2. Comparison of the nonlinear power 
spectrum in ACDM model with observational re- 
sults. Dots are results for the APM galaxy survey. 
Results for the real-space power spectrum for the 
CfA survey are shown as open circle (101/i _1 Mpc 
sample) and triangles (130/i _1 Mpc sample). For- 
mal error bars for each of the surveys are smaller 
than the difference between the results. The full 
curve represents the power spectrum of the dark 
matter. Lower limits on the power spectrum 
of galaxies predicted by the ACDM model are 
shown as the dashed curve (ACDM^ simulation 
in Ref. [5]) and the dot-dashed curve (ACDM C 
simulation). 
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In our recent paper [5] we show that lower lim- 
its on the predicted galaxy correlations in this 
ACDM model are also much higher than the 
Stromlo-APM and CfA2 data. It seems that 
none of the simplest and most attractive schemes 
for the distribution of galaxies in the model can 
give correlation functions that agree with obser- 
vations. Galaxies cannot follow the dark matter. 
Neither the power spectrum nor the correlation 
function allow this. The simplest biasing mod- 
els (halos with overdensity above 200 or density 
above any reasonable threshold) do not work ei- 
ther: the correlation function discrepancy on 2- 
3 Mpc scales can be reduced, but the situation on 
smaller scales gets even worse. 

These problems do not mean that all variants of 
the ACDM model are inconsistent with observa- 
tions, but it implies that the most attractive vari- 
ants with large age of the Universe, large Hubble 
constant, and relatively large cosmological con- 
stant are very difficult to reconcile simultaneously 
with the observed clustering of galaxies and with 
the number density of galaxy clusters. 

4. DISCUSSION 

Thus, although a h = 0.7, fi = 0.3 ACDM 
model normalized to COBE nicely fits the data 
on intermediate scales, including the abundance 
of rich clusters, the agreement of the linear power 
spectrum with the data on smaller scales leaves 
no room for the nonlinear effects that are surely 
important there. Our minimum N-body estimate 
of the power spectrum in this model for wavenum- 
bers k = (0.4 — 1) h~ x Mpc is at least a factor 
of 2 too high compared to the Cf A results and at 
least a factor of 3 too high compared to APM. 
This is true if we assume that the galaxies trace 
the dark matter without significant bias, but it 
is also true if we just assume that galaxy den- 
sity is any monotonically increasing function of 
dark matter density. It is true even in an extreme 
model in which we assume that this function is 
zero for regions of low density and constant when 
the density is above the threshold to produce even 
a small galaxy. The only way this model can be 
compatible with the data is for the galaxy forma- 
tion process to result in strongly scale-dependent 



anti-biasing on small scales, but this seems hardly 
compatible with the observations supporting hi- 
erarchical scaling and indicating that the number 
density of luminous galaxies is highest in regions 
where the dark matter is densest, such as clus- 
ters [5]. We suggest that ACDM models will have 
a better chance of agreeing with observations if 
they have higher Slo, lower h, and a tilted spec- 
trum of primordial fluctuations (n p < 0.9). 

It is instructive to compare the Slo = 0.3, 
h = 0.7 ACDM model that we have been consid- 
ering with standard CDM and with CHDM. At 
k = 0.5/i Mpc" 1 , Figs. 5 and 6 of Ref. [41] show 
that the Q, v = 0.3 CHDM spectrum and that of a 
biased CDM model with the same a% = 0.67 are 
both in good agreement with the values indicated 
for the power spectrum P(k) by the APM and 
CfA data, while the CDM spectrum with erg = 1 
is higher by about a factor of two. As Figure 3 
shows, CHDM with Q, v = 0.2 in two neutrino 
species [12] also gives nonlinear P(k) consistent 
with the APM data. 




k (hMpc -1 ) 

Figure 3. Comparison of APM galaxy power 
spectrum (triangles) with nonlinear cold parti- 
cle power spectrum from CHDM model consid- 
ered in this paper (upper solid curve). The dot- 
ted curves are linear theory; upper curves are for 
z = 0, lower curves correspond to the higher red- 
shift z = 9.9. 
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Aside from the Hubble constant, the main po- 
tential problem for CHDM appears to be forming 
enough structure at high redshift. Although as we 
mentioned above the indications that the amount 
of gas in damped Lyman a systems is starting 
to decrease at high redshift z > 3 seems to be 
in accord with the available data, the large ve- 
locity spread of the associated metal-line systems 
may indicate that these systems are more mas- 
sive than CHDM would predict (see e.g., [42]). 
Also, Steidel et al. [43] have found objects by 
their emitted light at redshifts z = 3 — 3.5 ap- 
parently with relatively high velocity dispersions 
(indicated by the equivalent widths of absorp- 
tion lines), which they tentatively identify as the 
progenitors of giant elliptical galaxies. Assuming 
that the indicated velocity dispersions are indeed 
gravitational velocities, Mo & Fukugita (MF) [44] 
have argued that the abundance of these objects 
is higher than expected for the COBE-normalized 
SI = 1 CDM-type models that can fit the low- 
redshift data, including CHDM, but in accord 
with predictions of the ACDM model considered 
here. (In more detail, the MF analysis disfavors 
CHDM with Slj, J> 0.2 in a single species of neu- 
trinos. They apparently would argue that this 
model is then in difficulty since it overproduces 
rich clusters — and if that problem were solved 
with a little tilt n p « 0.9, the resulting decrease 
in fluctuation power on small scales would not 
lead to formation of enough early objects. How- 
ever, if Slj, « 0.2 is shared between two species 
of neutrinos, the resulting model appears to be 
at least marginally consistent with both clusters 
and the Steidel objects even with the assumptions 
of MF. The ACDM model with h = 0.7 consis- 
tent with the most restrictive MF assumptions 
has fi > 0.5, hence t < 12 Gyr. ACDM mod- 
els having tilt and lower h, and therefore more 
consistent with the small-scale power constraint 
discussed above, may also be in trouble with the 
MF analysis.) But in addition to uncertainties 
about the actual velocity dispersion and physi- 
cal size of the Steidel objects, the conclusions of 
the MF analysis can also be significantly weak- 
ened if the gravitational velocities of the observed 
baryons are systematically higher than the gravi- 
tational velocities in the surrounding dark matter 



halos, as is actually likely to be the case at low 
redshift for large spiral galaxies [45], and even 
more so for elliptical galaxies which are largely 
self-gravitating stellar systems in their central 
regions. 

There is another sort of constraint from ob- 
served numbers of high-redshift protogalaxies 
which disfavors ACDM. The upper limit on the 
number of z > 4 objects in the Hubble Deep 
Field (which presumably correspond to smaller- 
mass galaxies than most of the Steidel objects) is 
far lower than the expectations in low-Slo mod- 
els, especially with a positive cosmological con- 
stant, because of the large volume at high red- 
shift in such cosmologies [46]. So evidence from 
high-redshift objects cuts both ways, and it is too 
early to tell whether high- or low-Slo models will 
ultimately be favored. But one sort of constraint 
on CHDM models is likely to follow both from the 
high-z data just discussed and from the prelimi- 
nary indications from the recent results on cosmic 
microwave anisotropies at and beyond the first 
acoustic peak [47]: viable CHDM models cannot 
have much tilt (i.e., n p ~ 1). 
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